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Abstract: We consider hydrodynamics of M = 4 supersymmetric SU(N C ) Yang-Mills 
plasma at a nonzero density of i?-charge. In the regime of large N c and large 't Hooft 
coupling the gravity dual description involves an asymptotically Anti- de Sitter five- 
dimensional charged black hole solution of Behrnd, Cvetic and Sabra. We compute 
the shear viscosity as a function of chemical potentials conjugated to the three U(l) C 
SO(6)r charges. The ratio of the shear viscosity to entropy density is independent of 
the chemical potentials and is equal to 1/Att. For a single charge black hole we also 
compute the thermal conductivity, and investigate the critical behavior of the transport 
coefficients near the boundary of thermodynamic stability. 
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1. Introduction 

Recent studies of strongly coupled thermal gauge theories in the framework of the 
gauge-gravity duality ( [1], [2], [3], for a review see [4]) suggest that in all those theories 
in the regime described by gravity duals the ratio of the shear viscosity to volume 
entropy density is universal and equal to l/An [5], [6], [7]. As this result was obtained 
assuming zero densities of conserved charges, a natural question to ask is what happens 
when the chemical potentials conjugated to these charges are turned on. 

The simplest ten-dimensional gravitational background corresponding to a non- 
zero chemical potential in a dual four- dimensional finite-temperature field theory is the 
one of spinning near-extremal three-branes [8], [9], [10]. The number of independent 
commuting angular momenta that can be given to the three-branes is equal to the 
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rank r = 3 of the isometry group 50(6) of the space transverse to the branes. Upon 
dimensional reduction on S 5 one obtains the background corresponding to the five- 
dimensional asymptotically AdS black hole of the Reissner-Nordstrom type with three 
U(l) charges proportional to the angular momenta of the branes [10]. This background 
was found by Behrnd, Cvetic and Sabra [11] as a particular solution to the equations of 
motion of the five- dimensional M = 2 gauged supergravity. The solution corresponding 
to the dimensional reduction of the spinning three-brane metric has a translationally 
invariant horizon. (It can also be regarded as a black hole with a spherical horizon in 
an infinite volume limit.) 

In the AdS / CFT correspondence the isometry group of S 5 is interpreted as the R- 
symmetry group of the dual M = 4 supersymmetric Yang-Mills (SYM) theory. Three 
independent chemical potentials /ij can be introduced as the Lagrange multipliers to 
the three U(l) charges in the Cartan subalgebra of SO(6)r. Thermodynamics of the 
R-charged black holes in the context of the AdS/CFT correspondence was originally 
studied in [8], [12], [9], [10]. One feature relevant for our discussion is the thermody- 
namic instability 1 occurring for black holes with excessively large charge (or equivalently 
for the three-branes rotating too fast) [8], [12]. The instability may signal the onset 
of a phase transition. However, both in gravity and in the dual strongly coupled field 
theory picture it is not clear what the new phase might be. 

In this paper we study the hydrodynamic regime of the four-dimensional Af = 4 
SU(N C ) SYM theory with three non-zero chemical potentials in the limit of large N c and 
large 't Hooft coupling g\ M N c . Using Kubo formula, we compute the shear viscosity as 
a function of the three charges (or chemical potentials) and show that for any values of 
the charges in the thermodynamic stability domain the ratio of the viscosity to entropy 
density is equal to 1/47T. For a technically simpler case of a single charge black hole 
we explicitly compute thermal correlation functions of the stress-energy tensors and 
R-currents in the shear channel in the hydrodynamic approximation. The correlators 
exhibit a diffusion pole with the dispersion relation that confirms the result for the shear 
viscosity found from the Kubo formula. We also compute thermal conductivity. The 
ratio of thermal conductivity and shear viscosity obeys a simple relation reminiscent 
of the Wiedemann-Franz law for the ratio of the thermal conductivity to the electrical 
conductivity. Finally, we investigate the behavior of the transport coefficients near the 
boundary of thermodynamic stability and compute the corresponding critical exponent. 

The paper is organized as follows. We review the STU-model solution of Behrnd, 
Cvetic and Sabra in Section £| The shear viscosity for the three-charge black hole 
solutions is computed in Section ^. Starting from Section |, we specialize to the single- 

1 The related instability of the gravitational background was studied in [13]. 
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charge background: in Section |] we compute the correlation functions of the shear 
mode components of the stress-energy tensors and R-currents in the hydrodynamic 
approximation, in Section |^ we obtain the thermal conductivity, and in Section |] we 
discuss the critical behavior of the transport coefficients. For convenience of the reader, 
in Appendix [A] we outline the rescaling procedure necessary to obtain an asymptotically 
AdS charged black hole with a translationally invariant horizon from the black hole 
with a spherical horizon. In Appendix P we review relativistic hydrodynamics at finite 
chemical potential. 



2. The R-charged black hole background 

The STU-model solution to equations of motion of D = 5 Af = 2 gauged supergrav- 
ity was found by Behrnd, Cvetic and Sabra [11]. The relevant part of the gauged 
supergravity effective Lagrangian is 

R + T2 V ~ iGijF'F^ - Gyd^PX* + -^=e^ x e ljk Fl v F^A\ , (2.1) 



where F % , i = 1,2,3 are the field-strength tensors of the three Abelian gauge fields, 
X 1 are three real scalar fields subject to the constraint A 1 A 2 A 3 = 1. The metric on 
the scalar manifold is given by 

G l3 = idiag [(A 1 )" 2 , (A 1 )" 2 , (A 1 )- 2 ] . 



The scalar potential is 



3 1 



i=l 



The three-charge non-extremal STU solution is specified by the following background 
values of the metric 

ds 2 = _^-2/3 fk df 2 + H l/3 ( f -l dr 2 + r 2^2J ? (2 2) 
2 

f k = k-^ + ^-H, Hi = l + ^ } H = H 1 H 2 H 3 , (2.3) 
as well as the scalar and the gauge fields 
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The parameter k determines the spatial curvature of dQ 2 k : k = 1 corresponds to 
the metric on the three-sphere of unit radius, k = - to the metric on R 3 . It was 
shown in [10] that the k = solution arises as the Kaluza-Klein reduction on S 5 of 
the ten-dimensional metric describing spinning near-extremal three-branes. The three 
R-charges qi are related to the three independent angular momenta in ten dimen- 
sions. Since, strictly speaking, hydrodynamic regime is meaningful only in the case of 
a translationally-invariant horizon, in this paper we set 2 k = and 

dtt 2 3fi -> (dx 2 + dy 2 + dz 2 ) . 

Introducing the new radial coordinate u = r+fr 2 , where r + is the largest root of the 
equation /(r) = 0, we write the background fields in the form 

ds 2 b = -n- 2 ^^^- f dt 2 + n 1/3 ^^- (dx 2 + d y 2 + dz 2 ) + n^^-Au 2 , (2.5) 

3 

f(u) =H(u)-u 2 T[(l + K i ), Hi = l + K iU} Ki = \, T = r + /7tL 2 . (2.6) 
i=i r + 

The scalar fields and the gauge fields are given by 3 

* = wn-v A ^ = iwrv ^ = #na + ^) 1/2 - w 

Hi(u) LHi{u) L aj. 

The Hawking temperature of the background (|2.5|) is given by 

2 V /(1 + K 1 )(1+K 2 )(1 + K 3 ) °' 

The volume density of the Bekenstein-Hawking entropy is 

i=l 

where V3 is the spatial volume along the three infinite dimensions of the horizon, G5 = 
7rL 3 /2N 2 . The energy density and pressure are given by 

e = _^LJ](1 + Ki ), (2.10) 

i=l 
3 



7l 2 N 2 Tft t-t 

P = ^-JLl[(l + Kl ). (2.11) 
' i=i 



2 The background with a translationally-invariant horizon and related thermodynamics can also be 
obtained by taking an infinite volume limit of the k = 1 solution (see Appendix [A|) . 

3 We change the normalization of the gauge fields by a factor of y2/L. This normalization is used 
in the rest of the paper. 
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The densities of physical charges are 

7T r\kiL ttN 2 T* 



Pi 



4G 5 y/2V 3 



«1 



,1/2 



(2.12) 



1=1 



The chemical potentials conjugated to pi are defined as 



(Mi = A*(u) 



7tT a/21^ 



/2 



(2.13) 



For the grand canonical ensemble, where the system in thermodynamic equilibrium 
is characterized by the values of temperature and chemical potentials regulating its 
interaction with the surrounding heat bath, the appropriate thermodynamic potential 
is the Gibbs potential Q, 



^/V 3 = -P = 8- T H S - ^ PiPi ■ 

1=1 

The first law of thermodynamics reads 

3 

dP = sdT H + Pidfii . 

i=l 

A stable thermodynamic equilibrium is determined by the conditions 

V /T,jUj fixed ' V 'T, in fixed 



(2.14) 



(2.15) 



(2.16) 



The stability condition ( 2.16|) translates into the equation 

d 2 e(s,pi 



det 



dsdpi 



> 0. 



(2.17) 



Since Ki = 8ii 2 p 2 /s 2 and 

e(s,p*) 



3s 4 / 3 



n i+ 



3 2 2\ 1/3 
^ Pi 



2(27TiV) 2 /3 11 

the condition of thermodynamic stability implies the following constraint on 

2 — K\ — K,2 — ^3 + &1 K>2 ^3 > . (2-18) 
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3. Shear viscosity 



The simplest way to compute shear viscosity from the dual gravity background is to 
use Kubo formula. Kubo formula relates the shear viscosity to the correlation function 
of the stress-energy tensor at zero spatial momentum, 

77 = lim -L / dtdxe™ ([T xy (x), T xy (0)}) = - lim ImG ^°) , (3.1) 

where the retarded Green's function for the components of the stress-energy tensor is 
defined as 

GW>,<7) = -i[ Ae-"^(t)([^(x), T Ap (0)]> . (3.2) 



Thus finding the shear viscosity amounts to computing the zero-frequency limit of the 
imaginary part of the retarded correlator G xyxy (u,q). To compute the correlator, we 
follow the procedure outlined in [14] and used in [15] to determine the shear viscosity of 
the strongly coupled Af = 4 SYM at zero chemical potential. First, one has to determine 
the equation obeyed by the component h xy (u, t, z) of the gravitational perturbation of 
the background ( |2.5| ) - (|2.7|) . By symmetry argument [15], [16] or by the direct analysis 
of perturbations of the equations of motion following from the Lagrangian ( |2.1| ) one 
can show that the Fourier component of the off-diagonal perturbation = h y decouples 
from all other perturbations and obeys the equation for a minimally coupled massless 
scalar in the background 



( t ) k + ^ ( t ) k + T2 ^fc = ' 3 - 3 

where 

W= 2^' q= 2^- (3 ' 4) 

The solution to Eq. (|3.3| ) in the hydrodynamic regime vo 1, q 1 can be obtained 
along the lines of Ref. [15]. We find 




iU_ 

T 



log 



c 3 (u 



a(u - l) 2 + b(u - 1) +c 

/h 2 — Anr\ 

+ ■■■ , (3.5) 



(2o(u - 1) + b - Vb 2 -4ac)(b + Vb 2 - 4ac) 



Vb 2 - 4ac (2a(u -l) + b + Vb 2 - 4ac)(6 - Vb 2 - 4ac) 
where ellipses denote higher order terms in the hydrodynamic expansion 



3 3 



J J Ki , b = 2a — 1 — /tj, c = b — a — 1 . 



i=l i=l 
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and 

£/(«i, « 2> «3) = f~ ; ; • 

In the limit Ki — > the solution ( |3~5l) reduces to the one found in [15]. 

Another essential ingredient in computing the correlator is the boundary action. 
The total action is given by 

S = — *— / d 5 xC + -^—[ d 4 xV^hK+-^— [ d A xV^hW, (3.6) 
1o7tG 5 J Mb 87tG 5 J dMs 87tG 5 J dMs 

where the second term is the Gibbons-Hawking boundary term, and the third term is 
required to make the action finite 4 in the limit u — > 0. The explicit form of W for the 
background of interest was determined in [21]: 

7/1/3 3 

w = -^ r Y, H i 1 - ( 3 - ? ) 

8=1 

Computing the action ( |3~6"1) on shell and expanding to quadratic order in h v x we find 
the following boundary action for <ftk 

5'b = -^2^^(ti)M«) • (3-8) 



Following the prescription of [14] the retarded correlator is then computed as 

GM= ^^mm. (3 . 9) 

(The solution <pk{u) has been normalized to 1 at u — e.) We find 

G xyxy {uj,q) = U{k 1 ,k 2 ,k 3 ). (3.10) 

o 

The Kubo formula ( |3.1| ) gives the shear viscosity 



3 



nN 2 T 3 H n(l + «i) 



\2 



77 = ^N 2 Tq v /(l + Kl )(l + K2 )(l + K3 ) = i=l . (3.11) 



(2 + «i + k 2 + k 3 - K1K2K3) 3 



4 For a discussion of relevant issues in the context of holographic renormalization, see e.g. [17], [18], 
[19], [20], [21]. 
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Comparing this result with the expression fl2.9|) for the entropy density we immediately 
conclude that for any value of the chemical potential 

2=1 (3.12) 

s 4tt v ' 



For small we have 



4. Shear viscosity from the diffusion pole 

In this Section we explicitly compute the retarded correlation functions of the stress- 
energy tensor and the R-currents and show that they exhibit a diffusion pole predicted 
by hydrodynamics. The value of the shear viscosity extracted from the pole is in 
agreement with Eq. (|3.11|) . For simplicity, in the rest of the paper we restrict ourselves 
to the case of a single charge black hole. We set q\ ^ 0, g 2 = q 3 = and omit the index 
"1" in all subsequent expressions. 

4.1 The single charge black hole background 

For a single charge black hole the effective Lagrangian (|2.1|) can be written as 



£ R+^V- ^-H^F 2 - iH^g^d.Hd.H, (4.1) 



L 2 8 3 

where F^ v is the field-strength tensor of a U(l) gauge field, and V is the potential for 
the scalar field H, 

V = 2H 2/3 + AH- 1/3 . (4.2) 

The system of the gauged supergravity equations of motion for the fields g^ u , A^, H 
reads 

UH = H-Y v d»Hd„H + ^H 7 / 3 F 2 - |^ 2 f^ ' 

d^yf^H^F^) =0, (4.3) 



2 



R/j-u = -^H^F^Fv 7 + ^H^ 2 d fl Hd u H — g^ u 



2 . . L 



2 



3L 2 24 

Consider small perturbations of the single charge background 



V+ t^ H ^F 2 



g^u = g% + h ^(u,t,z) , (4.4) 

A M = Al + A^u,t,z), (4.5) 
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where g® v , A^ are given by Eqs. (|2.5| ), ( p.7|) with K\ = k, k 2 = k 3 = 0. We assume 
that the perturbations depend on time, radial coordinate, and only one of the spatial 
world-volume coordinates, z. We are interested in gravitational fluctuations of the 
shear type, where the only nonzero components of h^ u are h ta , h za , a — x,y. One can 
show that fluctuations of all other fields except A a (r, t, z), a — x, y, can be consistently 
set to zero. Introducing the new variables 

2A 

T = H ta = g xx h ta , Z = H za = g xx h za , A=^ (4.6) 
the system of linearized equations derived from Eqs. (|4.3|) can be written as 



uH ju ju 2H 

„ uf — f , w 2 H o w H , 
Z" + J - J Z' + — - Z + 5— T = , 4.7c 

uf f 2 u f 2 u 

a 2 H va 2 H 2 

(HfA' + 2(1 + k)T)' - - — A + ^— A = 0. (4.7d) 

u ju 

These equations are not independent: combining Eq. (|4.7a|) with Eq. ( |4.7b|) , one obtains 
Eq. (|4.7c|) . Thus it is sufficient to consider Eqs. ( |4.7a|) , ( |4.7b|) , (|4.7d|) . Expressing 
Z(u) from Eq. ( [4.7b| ), we differentiate it with respect to the radial coordinate u and 
substitute the resulting expression for Z'(u) into Eq. ( |4.7a| ). Thus we obtain a system 
of two coupled differential equations for G{u) = T'(u) and A(u) 

H ' n c , ( H ' f' ™ 2 h q 2 



G " + \H Jr l) G ' + \^H 2 ~^fH + ^P~-f) G 

nu KuHf> + Kf(2H-uH>) + ^ A = Q (4g) 
2H 2fH 2 2f 2 ' V ; 

(For k = Eqs. (fl8l) , ( [4.9|) decouple 5 and reduce respectively to Eqs. (6.15) and (5.5d) 
of Ref. [15].) For the system ( |4.8|) , (|4.9|) , the exponents at the singular point u = 1 
corresponding to the horizon are 

a = { , , ivoU(k) , -ivoU(k) } , (4.10) 



3 One should recall the rescaling (|4.6| 
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where U(k) = U(n,0,0). Choosing the exponent —ixx>U(K.) corresponding to the in- 
coming wave boundary condition at u = 1, we look for the solutions to Eqs. (|4.8|) , ( |4.9|) 
in the hydrodynamic approximation in the form 

G(u) = - U{ * + K $ r mU{K) (G (u) + xoG^u) + q 2 G 2 (n) + •■■), (4.11) 



where functions Gi, A{ are regular at u = 1. We obtain 

G (u) =d-—, 

K 

A (u) = d + uC 2 , 
where C\, C 2 are the integration constants. Next, 

i 



K 1 (k + 2) u (2 + ku) 

— G\ K (2 + K + KU + K 2 u) 



w-1) C 2 (k + 2) (2 + ku) 



(4.12) 

(4.13) 
(4.14) 



+ 2 (k + 1) (C 2 + uCi) it (2 + ku) log 



k + 2 



0c + 1 (k + 2) 
+ 2 (« + 1) (d + C 2 m) log 



u + #(«) j ' 
« - 1) ^Ci k - C 2 (« + 2) 

K + 2 



u + H(u) 



(4.15) 



(4.16) 



Functions G 2 (u), A 2 (u) are given by rather cumbersome expressions involving polylog- 
arithms. We do not write them explicitly here. 

The integration constants C±, C 2 can be expressed in terms of the boundary values 
of the fields , , by solving the equations 

UmA{u) = A (0) 

1 



u->0 
KU 



li mu f(u)(G'- — G , , 



A') =t»qZ(°»+q 2 T(°) 



(4.17) 
(4.18) 



where Eq. ( 4.18Q comes from Eq. ( 4.7b| ). We find 

C x = kA<®, 

c 2 -- 



2k(\ + k)< [ 2 t (q) K 2 (q 2 -iVl + Kn>) ^ (0) 



q 2 — i 2 vT + re to 



q 2 — i 2 + «; to 



(4.19) 
(4.20) 
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One can observe the appearance of the hydro dynamic pole in Eq. (|4.20 ). 



4.2 The correlators 

We are interested in two-point retarded correlation functions of stress-energy tensors 
and i?-currents defined by 



G^ Xp {u,q) = -i \ d 4 xe-^(t)([7V(x), T Ap (0)]> . 



(4.21) 



(4.22) 



(4.23) 



To compute the correlators we need to consider the boundary action. On shell, the 
action reduces to the surface terms, S = Shorizon + S t , where 



G^xfaq) = -i / cfxe-^(t)([T^(x), J A (0)]> . 



G^(u,q) = -z/^e-"^(t)([J^), J„(0)]) 



u^o 8 



-(1 +k)- ~H' ta H ta + -H' za H za + ^- 2 A' a A a 



II 



U 



~ f (1 + k) K - 1(1 + «) Hi + V2 f^ K) H ta A a 



(4.24) 



The first term in Eq. fllgD , -7r 2 iV 2 T 4 (l + k) V4/8, is the density of the Gibbs potential 
f2 (i.e. the pressure with a minus sign), times the four- volume. The retarded two-point 
functions are obtained from S e following the recipe formulated in [14] . After substitut- 



ing the solution ( |4.12|) into Eq. (|4.24D , the part of the boundary action quadratic in 
fluctuations assumes the form 



dudq 



(4.25) 



where 0- ' 1 denote the boundary values of the fields T^°\ ■ Then the retarded 

correlators are given by 6 



G 1 



-2F ik {u,q), i = k, 



(4.26) 



6 Note that we obtain the correlators with upper indices. Indices of the boundary theory correlators 
are raised or lowered with the flat Minkowski metric, so that e.g. G taza — —Gtaza- 
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Computing the correlators, to leading order in the hydrodynamic approximation we 
obtain 7 

Gtata(uJ, q) = ——. — ^- , 4.27 

— Vq 2 ) 



VTT^iV 2 7rT 3 o;g 
8(icu - Vq 1 



G taza {u, q) = - Q/ , . - ° , (4.28) 



G zaza {u, q) = — — , 4.29 

where the diffusion constant is given by 

1 1 I + k/2 , 

AnToVTT^ 4nT H 1 + k y 1 

In the limit k — > the results (|4.27|) - (|4.30Q reduce to those obtained in [15] for the case 
of a zero chemical potential. For the correlators of the components of a stress-energy 
tensor and an .R-current we have 



Gtaa{V, q) = ,. • 4.31 

6{tio — Dq A ) 

„ . N V2^N 2 T 2 uq 

G zaa (", 9) = ~^f. ±4 • 4.32 

62{iuj — Vq 2 ) 

These correlators vanish in the limit k — > 0, in agreement with [15]. Finally, the 
retarded correlator of the x (or component of the R-currents is given by 

ikN 2 T 2 uj 

G xx (lu, q) = G yy (u, q) = G aa (u, q) = — °— + 0(u, q 2 ) . (4.33) 

IQ(iuj — Vq 1 ) 



In the limit k — >• the leading contribution in Eq. ( |4 . 3 3| ) vanishes. The subleading 
term gives G aa = —iN 2 Tquo / IQn which coincides with the result obtained in [15]. 

In the limit of vanishing spatial momentum the nontrivial contribution to G aa (u), q) 
again comes from the subleading term in Eq. ( |4.33j ) . It is given by 

„ , \ i(K + 2) 2 N 2 T LU 

G " W ' 0) = ~ RA I — T"T • ( 4 - 34 ) 
647TVK + 1 

We also find that in the limit q — > the correlators Gtaa and Gtata vanish (modulo 
contact terms). 



7 Contact terms are ignored. 
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4.3 The diffusion pole 

All the retarded correlators in the shear channel exhibit a diffusion pole with the 
dispersion relation 

uj = -iVq 2 , (4.35) 

where the diffusion constant T> is given by Eq. ( }4.30|) . To find viscosity, recall that in 
hydrodynamics 

v 77 



e + P 

From thermodynamics it follows that 

2(1 + k) m 
e + P = T H s + pp = -4- >- T H s . 

2 + K 

Thus for the ratio of shear viscosity to entropy density we find rj/s = 1/4tt which 
coincides with the result ( |3.12| ) obtained from the Kubo formula. 



5. Thermal conductivity 

Thermal conductivity kj- can be computed using the appropriate Kubo formula (see 
Appendix |B]) 

(e + P) 2 1 

k t = - V > lim - Imffw, . (5.1) 

Here G is the retarded Green's function of the R-current components J x given by 
Eq. flPD . Thus we find 

32tt p 2 h k(k + 2) v ; 

In terms of the chemical potential m the thermal conductivity can be written as 



KT = ^-r H i + v / ^?-^(v^w-5) 

8m 

Comparing this result to the one for the shear viscosity ( |3.11| ), we observe that for a 
single-charge black hole, the shear viscosity and thermal conductivity can be expressed 
in terms of the chemical potential pi = p as 

TtN 2 T% , s -nN 2 T 2 . . . , 

7i=—^ JL F ri (ji,T H ), K T = —-^F K (p,T H ), (5.4) 
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where the functions F v , F K depend only on the ratio m = h/27tTh, 

8m 2 (l - v 7 ! -4m 2 -m 2 ) 2 2 

F V (H,T H ) = \ -3 J , F K (fi,T H ) = —F v (fi,T H ). (5.5) 

(l - v / T T W) m 

Thus for all values of T H and \x one finds an analogue of the Wiedemann- Franz law [24] 



2 



For small m we get 



= 8tt 2 . (5.6) 



F v = l + m 2 -m 6 + 0(m 8 ). 
The function F v (m) is shown in Fig. 1. 

6. Critical behavior of transport coefficients 

The boundary of thermodynamic stability is m c = 1/2 or fi c = ttT. Expanding the 
shear viscosity and the thermal conductivity near m = m c we obtain 



V = V* 
kt = kt 

where 



2 20 
1 + - vV-m - — (m c - m) + O ((m c - m) 3/2 ) 

16 



1 + - ^m c -m + — (m c - m) + O ((m c - m) 3/2 ) 



y7riV 2 T 3 9nN 2 T* 

V* — 7T A , «T* " 



(6.1) 
(6.2) 



64 1 8 
Both the viscosity and the thermal conductivity are finite 8 at the critical point. Their 
derivatives diverge with the critical index equal to 1/2. 

7. Conclusion 

We have considered the hydrodynamic regime of the M = 4 supersymmetric theory at 
finite temperature and finite chemical potential. We have computed the shear viscosity 
and the thermal conductivity. The shear viscosity is computed using two different 
methods which give the same answer. We find that the ratio of shear viscosity and 
the entropy density is always equal to 1/47T, which is the same value found so far in 



8 Curiously, the shear viscosity of He 4 near the A-point is also finite and its first derivative is 
divergent, as first shown theoretically by A. M. Polyakov [25] and later confirmed experimentally [26]. 
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Figure 1: Normalized shear viscosity F n as a function of the chemical potential m = h/2-kTh- 
There is a cusp singularity at \i c = ttTh- The part of the curve to the right of the singularity 
(shown in dashed line) is unphysical. 

all theories with gravity duals. Our result demonstrates that the universality of this 
ratio extends to theories with gravity duals at finite chemical potentials. We found a 
curious relationship between the shear viscosity and the thermal conductivity similar 
to the Wiedemann-Franz law, and we have also determined the critical behavior of the 
kinetic coefficients near the boundary of thermodynamic stability. 

One possible extension of this work is to compute the thermal conductivity and dif- 
fusion coefficients in the case when all three conserved charges are nonzero. This would 
complete our knowledge of the kinetic coefficents of M = 4 super- Yang-Mills theory 
in the whole phase diagram, as the conformal invariance of the theory guarantees that 
the speed of sound is equal to 1/ v^3 and the bulk viscosity is zero for any temperature 
and chemical potential. 

Note added: while this paper was being completed, we became aware of the work 
on the same subject by J. Mas [27] whose results are in agreement with our analysis. 
Another recent paper on the subject is by K. Maeda, M. Natsuume, and T. Okamura 
[28]. A closely related work on the hydrodynamics of M2-branes by O. Saremi appeared 
in [29]. 
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A. Rescaling the black hole solution 

The metric (2.5) of the gravity background dual to Af = 4 SYM with nonzero chemical 



potentials in Minkowski space can be obtained by rescaling and taking an infinite vol- 
ume limit of the black hole solution with a spherical horizon. (At zero chemical poten- 
tial, this procedure leads from the metric corresponding to an AdS-Schwarzschild black 
hole to the metric describing the near-horizon region of non-extremal three-branes.) 
Following Refs. [22] , [9] , here we explain how the rescaling works for the solution itself 
as well as for various thermodynamic quantities associated with it. 

Focusing for simplicity on a single charge black hole, the metric (|2.2|) with k — 1 
reads 

dsl = -H- 2 /'hdt 2 + H 1 '* (f^dr 2 + r 2 dQ 2 3 ) , (A.l) 



where 



Rescaling 



iy £i j_J ty^ 



r^A 1/4 r, t^\- 1/4 t, mi ^ Ami, q -> A 1/2 q , (A.2) 
and taking A —>■ oo while simultaneously blowing up the sphere 

L 2 dn 2 3 -> \- 1/2 (dx 2 + dy 2 + dz 2 ) (A.3) 
in the limit we obtain the metric with a translationally invariant horizon 

ds 2 5 = -H- 2/ ^ fdt 2 + H 1 ** ^ dr 2 + H^%- (dx 2 + dy 2 + dz 2 ) , (A.4) 

where / = 1 + q/r 2 — r 4 /r 4 , with r 4 = m\L 2 . 

Similar reasoning applies to thermodynamic quantities and their densities. Tem- 
perature, entropy, energy and the Gibbs potential scale as 

T H ^\ 1/4 T H , S^S, E^\ 1/4 E, n^\ 1/4 tt. (A.5) 
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For the background ( |A.1|) , the inverse Hawking temperature, entropy, energy, and the 
thermodynamic potential Q are given correspondingly by (see e.g. [21]) 

a o Tl ( r l + l) 1/2 a 2 / 2 , \l/2 

P = 27TL l + q + 2ry S BH = ^-r + (r + + q ) , 

e = — —^ + 1 + — ) , n = — -_JL + -± + 



G 5 V8L 2 4 32 J ' G 5 V 8L 2 4 32 

In the limit A — > oo for the rescaled quantities we find 

„ o r2 ( r + + g) V2 c X 3/4 71-2 2 / 2 , \l/2 

From (|A.3|) one can see that 2tt 2 L 3 — > A~ 3 / 4 V3 and thus the pressure and the densities 
of entropy and energy are finite in the A —>■ oo limit and are given by 

CA.U 

^=-^ = 16^- (A ' 9) 



With identifications 9 and G 5 = nL 3 /2N 2 from Eqs. flOD , (|A~8D , (|OD we obtain 
the expressions (|2.9j) , ( p.lOj ), ( [2.1 1|) used in the main text. 



B. Relativistic hydrodynamics at finite chemical potential 

For completeness here we review the hydrodynamics of a relativistic fluid with one 
conserved charge. The hydrodynamic equations include the continuity equations 

d^ v = 0, <V" = (B.l) 

and the constitutive equations, which formally have the form 

r=(e + P)!iV + Ff + r F , J» = pu fl + ^ (B.2) 



Note that 7q = ri_(ri_ + q), where r + is the largest root of the equation f(r) = 0. 
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Here e and P are the local energy density and pressure, is the local velocity, u^u 11 = 
— 1. The parts r^ u and are the dissipative parts of the stress-energy tensor and the 
current. To complete the system of equations we need expressions relating t^ v and v 11 
with derivatives of u M and of the thermodynamic potentials. 

Following Landau and Lifshitz [23], we can choose and p so that t^ v and are 
orthogonal to u M 

u^ v = = . (B.3) 

The most general form of the constitutive equation follows from the second law of 
thermodynamics. First we notice that 

u v d^ v = -(e + P)c>X - + u v d^ v = (B.4) 

Using the thermodynamic relations 

e + P = Ts + pp, de = Tds + pdp, (B.5) 

current conservation, and Eq. ( P-3|) , Eq. ( |B.4|) can be transformed into 

d tl {sv?) = ^d tl v»-—d tl u v (B.6) 

3, - ^) = -v»dj^ - T —d^ . (B.7) 

We now interpret the left hand side as the divergence of the entropy current. The right 
hand side thus must be positive. This implies 



■*(^ + ^)> (B.8) 



= -rj^u" + d u u" + u»u x d x u u + u v u x d x u») - ( ( - \ ) (g» v + uV)d x u x 



(B.9) 

Here t] and ( are the shear and bulk viscosities, respectively. To have an interpretation 
of x as the coefficient of thermal conductivity, let us consider the case when there is no 
charge transport, J 1 = 0, but there is an energy flow, T tl ^ 0, which is the heat flow. 
The local velocity u % is necessarily small and is equal to 

u *=*#£. (B.10) 
p T y 1 



Therefore 



T u = (e + PK = (e + P). (B.li; 
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Using dP = sdT + pdfi one can write this equation as 

T " = -*(!-w-)\ d ' T -7h d ' p )- (B12) 

In the nonrelativistic theory the heat flow is proportional to the gradient of temperature; 
in the relativistic limit there is an extra contribution proportional to the gradient of 
pressure. The proportionality coefficient is the thermal conductivity, 

« r =^£±^Vx. (B.13) 

The Kubo's formula for k t can be written down as follows. Suppose one puts 
the thermal system in a slowly-varying external background gauge field coupled 
to the conserved charge. This field will induce a current, proportional to the electric 
field Ei = d t Ai — d^A t . But since A t plays the same role as the chemical potential, by 
comparing with Eq. ( |B.8| ) we can write 

J* = ^(d t A i - <9 l 'A*) (B.14) 

In the case when the external fields are homogeneous in space, this relation becomes 
very simple, 

J l = i^uA\ (B.15) 

We can compare this relation with the one that follows from the linear response theory, 
J* = -G R (u, 0)A\ Thus we obtain 

G R (u,0) = -iu^. (B.16) 



From Eq. (B.13) we find the Kubo formula (|5.1|). 
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